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Abstract 

The inelastic interaction between heavy ions and an electron plasma in the presence of an intense radiation 
field (RF) is investigated. The stopping power of the test ion averaged with a period of the RF has been calculated 
assuming that cjo > ui p , where ujo is the frequency of the RF and oj p is the plasma frequency. In order to 
highlight the effect of the radiation field we present a comparison of our analytical and numerical results obtained 
for nonzero RF with those for vanishing RF. It has been shown that the RF may strongly reduce the mean energy 
loss for slow ions while increasing it at high- velocities. Moreover, it has been shown, that acceleration of the 
projectile ion due to the RF is expected at high-velocities and in the high-intensity limit of the RF, when the 
quiver velocity of the plasma electrons exceeds the ion velocity. 
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1 Introduction 

The interaction of charged particles with a plasma in the presence of radiation field ( RF) has been a subject of 
great activity, starting wi t h the work of Tavdgiridze, Aliev, Gorbunov and other authors ( Tavdgirid ze & Tsintsadze 



^dg 

1970; Aliev et all Il97lt lArista et all Il989t IXkopvan et all Il997t iNersisvan & Akopvanl Q999). A comprehen 



sive treatment of the quantities related to inelastic particle-solid and particle-plasma interactions, like scattering 
rates and differential and total mean free paths and energy losses, can be formulated in terms of the dielectric 
response function o btained from the electron gas model. The results have importa nt applications in radiation and 
solid-state physics ([Ritchie et all Il975t iTung & Ritchiei 1 19771 lEcheniquel 1 19871). and more recently , in studies 



of energy deposition by ion beams in inertial conf i nement fusion (ICF) targets (Arista & Brandt, 1981; Me hlhorn 
19811: Mavnard '&Deutschl 1 19821 lArista & Pirizi 1 19871: b'Avanzo et all ll993UCouiilaud et all Il994h . On the 



other hand, the achievement of high-intensity laser beams with frequencies ranging between the infrared and 
vacuum-ultraviolet region has given rise to the possibility of new studies of interaction pro cesses, such as electron- 
atom scattering in las er fields (Ki-oll & Watsonl 1 1 9731; fWeingartshofer et all 1 19771 1 1983b . multip hoton ionization 
jLompre et all Il976t iBal dwin & Bore hami 1 198 lb . inverse bremsst rahlung and plasma heating (ISeelv & Harris , 
1973; Kim&Pac, 1979: iLimaef oil Il979h . screening breakdown (Miranda et al. . l2005h . and other processes of 
interest for applicati ons in optics, solid-state, and fusion research. In addition, a promising ICF scheme has been 



recently proposed (Stock fef q/lll996t iRoth et al . 2001), in which the plasma target is irradiated simultaneously 



by intense laser and ion beams. Within this scheme several experiments ( Frank et al. , l2010tlHoffmann et a/.Ll2010b 



have been performed to investigate the interactions of heavy ion and laser beams with plasma targets. An important 
aspect of these experiments is the energy loss measurements for the ions in a wide-range of plasma parameters. It 
is expected in such experiments that the ion propagation would be essentially affected by the parametric excitation 
of the plasma target by means of laser i rradiation. This effect has been supported recently by particle-in-cell (PIC) 
numerical simulations jHu et a/.luOl II) . 

In this paper we present a study of the effects of intense RF on the interaction of nonrelativistic projectile 
ions with an electron plasma. Our objective is to study two regimes of the ion energy loss which have not been 
considered in detail. For the first part of our study, we consider energy loss of a slow ion. In particular, this is 
motivated by the fact that the alpha-particles resulting from the nuclear fusion in a very dense plasma with tem- 
perature in the keV range, display a velocity mostly below electron thermal velocity. The second objective of our 
study is to investigate the energy loss in high-velocity regime. Previously this has been done for a classical plasma 
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(Tavdgiridze & Tsintsadze, 1 19701: Aliev et al\ . I1971I: In ersisvan & Akopvarl 1 19991 ) treating only the collective ex- 
citations as well as in the rang e of solid-state densities (fully degenerate plasma) and at the intermediate intensities 
of the RF ( Arista et a/.lll989l) when the electron quiver amplitude is comparable to the screening length of the tar- 
get. To gain more insight into the RF effect on the energy loss process we consider here the regime of intense RF 
when the quiver amplitude largely exceeding the typical screening length of the fully degenerate electron plasma. 

The plan of the paper is as follows. In Sec. [2] we briefly outline the RPA formulation for the energy loss of a 
heavy ion uniformly moving in a plasma in the presence of an intense RF. The limiting case of a weak RF is also 
considered. In Sees. [3] and |4] we have calculated the effects of the RF on the mean energy loss (stopping power) 
of the test ion considering two somewhat distinct cases with slow (Sec. [3} and fast (Sec. H} projectiles moving 
in a classical and fully degenerated electron gas, respectively. In the latter case the d egenerated electron gas i s 
treated within a simple plasmon-pole approximation proposed by Basbas and Ritchie (IBasbas & Ritchie . 1982 ). 
It has been shown, that besides usual stopping in a plasma it is possible to accelerate the charged particles beam 
through RF. This effect is expected for fast projectiles and in the high-intensity limit of the RF, when the "quiver 
velocity" of the plasma electrons exceeds the projectile ion velocity. The results are summarized in Sec. |5]which 
also includes discussion and outlook. 



2 RPA formulation 

The whole interaction process of the projectile ion with a plasma involves the energy loss and the charge states 
of the ion and - as an additional aspect - the ionization and recombination of the ion driven by the RF and the 
collisions with the plasma particles. A complete description of the interaction of the ion requires a simultaneous 
treatment of all these effects including, in particular, the effect of the ion charge equilibration on the energy loss 
process. In this paper we do not discuss the charge state evolution of the projectiles under study, but concentrate on 
the RF effects on the energy loss process assuming an equilibrium charge state of the ion with an effective charge 
Ze. This is motivated by the fact that the charge equilibration occurs in time scales which are usually much smaller 
than the time of passage of the ion through target. 

The problem is formulated using the random-phase approximation (RPA), and includes the effects of the RF 
in a self-consistent way. The electromagnetic field is treated in the long-wavelength limit, and the electrons are 
considered nonrelativistic. These are good approximations provided that (1) the wavelength of the RF (Ao = 
2itc/(jJo) is much larger than the typical screening length (A s = v s /lu p with v s the mean velocity of the electrons 
and uj p the plasma frequency), and (2) the "quiver velocity" of the electrons in the RF (ve = eEo/muJo) is much 
smaller than the speed of light c. These conditions can be alternatively written as (1) i*jq/lu p <C 2ttc/v s , (2) Wl <§C 
^noc(mc 2 )(uJo / 'w p ) 2 , where Wl = cEq/8tt is the RF intensity. As an estimate in the case of dense gaseous 
plasma, with electron density no = 10 18 cm -3 , we get i«oTOc 3 ~ 1.2 x 10 15 W/cm 2 . Thus the limits (1) and (2) 
are well above the values obtained with currently available high-power RF sources, and so the approximations are 
well justified. 

We consider the time-dependent Hamiltonian for the plasma electrons in the presence of both a r adiation field 
(RF) with vector potential Aft) = (c/^q)Eq cos(woi), and a self-consistent scalar potential (p(r, t) dArista et al. 
19891; iNersisvan & Akopva"ni ll999). i.e., 



ff W =E^(P- e -A(t)) 2 c+c p - e5>(M) C + +kCp , (1) 

P P,k 

where c p , c+ are annihilation and creation operators for electrons with momentum p, respectively, and </?(k, t) is 
the Fourier transform of (p(r, t). 

The potential ^(k, t) is produced by the external charge and by the induced electronic density, viz., 

fcV(k, t) = 4 7 r jO0 (k, t) - Aire £ iV p (k, t) (2) 

p 

being po(k,t) the Fourier transform of the external charge density po(r,t), and N p (k,t) = (c p _ k c p )t is the 
electrons number operator. 

The time evolution of the operator Ap(k, t) is determined by the equation 



m dN p (k,t) =[JVp(M)>g(t)]< (3) 



In particular, for an oscillatory field A(t) and within random-phase approximation Eq. ((3) has the solution 
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dArista et al.Lll989tlNersisvan & Akopvanill999h 



dt'ip(k, t') exp 
x exp [— i£ (sin(wot) — sin(woi'))] i 



JVp(k,t) = -(/ p _ k -/ p 



^ (e P -k - E P ) (* - 



(4) 



where £ — k ■ a, a = eEo / mw§ is the oscillation amplitude of the electrons driven by the RF (quiver amplitude), 
e p = p 2 /2m is the electron energy with momentum p. Here f p is the equilibrium distribution function for the 
electron plasma. 

Finally, using Eq. (f2]i and making a further Fourier transformation we obtain a solution for the potential ip in 
the form 

47rp (k,o;) 



where we have introduced the frequency transforms tp(k, lu), po(k, uS) of the quantities 



A)(M) 



A)(M) 



i£ sin(w t) 



(5) 



(6) 



and e(fc, w) is the RPA dielectric function jLindhardlll954l:lLindhard & Wintheiill964 . 

We consider a heavy point-like particle with mass M and effective charge Ze which moves with rectilinear 
trajectory with constant velocity v. We thus neglect the effect of the RF on the particle assuming that the quiver 
velocity of the ion in the laser field v q = ZeEo /Mojq <C V s , v. Here v s is the mean velocity of the target 
electrons. The charge density of the point-like ion is then given by po(r, t) = Ze5(r — vt). Inserting the Fourier 
transformation of this formula with respect to r into Eq. © and making a further Fourier transformation we obtain 



Pol 



(k, uj) — 2irZe ^ J n (()fi — k • v + nuj ) 



(7) 



where J n is the Bessel function of nth order. Using Eqs. (H)-© for the self-consistent potential ip(r, t) we finally 
arrive at 



¥>(*.*)= 2^ 



J(n— rn)u} t 



m.n— — oo 



rfk 



e ik-(,-vt) J m (C) J n (C) 

k 2 e(k, k • v — nwo) 



(8) 



This result represents the dynamical response of the medium to the motion of the test particle in the presence of 
the RF; it takes the form of an expansion over all the harmonics of the field frequency, with coefficients J n (£) that 
depend on the intensity Wl oc a 2 . 

From Eq. © it is straightforward to calculate the electric field E(r, t) = — V(^(r, t), and the time average 
(with respect to the period 2tt/luq of the laser field) of the stopping field E stop = (E(vi, t)) acting on the particle. 
Then, the averaged stopping power (SP) of the test particle becomes 



c- r, v _ _ 2Z 2 e 2 
o = -Ze— • E s top — /„ x 2 / 



oo 

E 

oo 



dk kv ^ 2 



k 2 



J 2 (C)Im. 



■1 



(9) 



with n n (k) = nu)Q + k • v. 

To illustrate the effects of the RF it is convenient to take into account the symmetry of the integrand in Eq. (0, 
with respect to the change k, n — > — k, —n. Using also the property of Bessel functions, J 2 n (C) = >^n(C)> we 
obtain 



2n 2 v 



dk- 



k- v 



J 2 (C)Im 



-1 



s(k, k ■ v) 



2J2 ^(C)Im 



-1 



e(k,n n (k)) 



(10) 



Hence, the SP depends on the particle velocity v, the frequency uio and the intensity Wl = cE 2 / 8ir of the RF 
(the intensity dependence is given through the quiver amplitude a). Moreover, since the vector k in Eq. (TTOb 
is spherically integrated, S becomes also a function of the angle between the velocity v, and the direction of 
polarization of RF, represented by a. 

By comparison, the SP in the absence of the RF is given by ( Deutschi 1986t Peter & Mever-ter-VehnI 1991 ) 



Sr = 



Z 2 e 2 
2n 2 v 



dk——— Im 



-1 



e (k, k ■ v) 



(11) 
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In the presence of the RF the SP Sb is modified and is given by the first term in Eq. (TlOb ("no photon" SP) 



S = 



2-k 2 v 



k 2 



^o(C)Im 



-1 



e (fe, k • v) 



(12) 



Next we consider the case of a weak radiation field (a < A s , where A s is the characteristic screening length) at 
arbitrary angle § between v and Eo. In Eq. ( TTOb we keep only the quadratic terms with respect to the quantity a 
and for the stopping power S we obtain 



Z 2 e 2 

47T 2 l 



S — Sb 



dk 



(k-v)(k-a) 2 Im 



1 



1 



e(fc,wo + k-v) e{k, k-v) 



(13) 



where Sb is the field-free SP given by Eq. (fTTT >. Note that due to the isotropy of the dielectric function e(k, lu) the 
angular integrations in Eqs. (TlOb — (TT3T> can be easily done. 

It is well known that within classical description an upper cutoff parameter fc max = l/r m ; n (where r m ; n 
is the effective minimum impact parameter) must be introduced in Eqs. ( fTTT i and ( fT3l ) to avoid the logarithmic 
divergence at large k. This divergence corresponds to the incapability of the classical perturbation theory to 
treat close encounters between the projectile particle and the plasma electrons properly. For r m j n we use the 
effective minimum impact parameter excluding hard Coulomb collisions with a scattering angle larger than 7r/2. 
The resulting cutoff parameter fc max ~ m(v 2 + v 2 u )/\Z\e 2 is well known for energy loss calculations (see, e.g., 
Zwicknagel et al. dl999h : Nersi svan et al.\ (2007) and references therein). Here w t h is the thermal velocity of the 
electrons. In particular, at low projectile velocities this cutoff parameter reads fe max = T/\Z\e 2 , where T is the 
plasma temperature given in energy units. 



3 Energy loss of slow ions 



In this section subsequent derivations are performed fo r the classical plasma and in the low-velocity limit of the 
ion. In this case the RPA dielectric function is given by (iFried & ConteL ll961) 



e(k,oj) = 1 + tJt^-W ' 



fc 2 A D \ fcl, th / ' 



(14) 



where Ad is the Debye screening length, and W{z) — g(z)+if(z) is the plasma dispersion function (IFried & Conte 
1961b with 



g(z) = 1 - ze 



(15) 



Consider now the SP determined by Eq. (ITOb in the limit of low- velocities, when v -C Wth- As discussed above 
we also assume that v 3> v q and neglect the effect of the RF on the ion. In the limit of the low-velocities from 
Eqs. (ITOb— (IT5l> we obtain 

S(7,a,0) = £1,3(7,0,0), (16) 



where 



5(7, a, ■&) = Ex (7, a) + S 2 (7, a) sin 2 0, 



(17) 



S a (7,a) 



e k 3 dk 



m Jo (k 2 + l) 2 Jo 



Jo(AkfJ,)f s (fJ,)dfi 



(18) 



Wig/ - 



Wxin/k^dk 



(k 2 +W(n/hy)) s 



JliAk^UliW] 



Here s = 1, 2, and = /i 2 , /2(m) = |(1 — 3/i 2 ). Note that at the absence of the laser field (i.e., at a — > 0) 

Si(7,a) — » 1, 52(7,0) — » 0. In this case the SP is determined by the quantity Sb in Eq. (TTTb dDeutsch , 1986t 



Peter & Mever-ter- Vehnl Tl 99 lb 



2Z 2 e 2 v 

7T 6X1. V th 



where 



(19) 



(20) 
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fl/A D 

Figure 1: The dimensionless quantities H(7, a, i9) (the lines with symbols) and E av (7, a) (the solid line without symbols) vs 
the intensity parameter of the laser field a/ Ad for ■£) — (solid line), ■& — 7r/4 (dashed line), $ = n/2 (dotted line) and for 
LJQ = 1.2o?p, 



is the Coulomb logarithm with £ = fc max AD- Also in Eqs. (TT6b — (fl~8T> we have introduced the angle -d between the 
velocity v and the polarization a vectors, W\(z) = dW(z)/dz, A = cl/\d, 7 = w p /wo < 1. Note that while the 
k integral in Eq. (fTTT i diverges logarithmically in a field-free case, Eqs. (fTZl i and ( f]~8T > are finite and do not require 
any cutoff. The Bessel functions involved in these expressions due to the radiation field guarantee the convergence 
of the fc-integrations. However, since in the sequel we shall compare Eqs. dT6b— (TT8l> with field-free SP Sb, for 
consistency the upper limits of the /c-integrals in Eq. (\M are kept finite with the same upper cutoff parameter as 
in Eqs. <fTTT > and ( Tl9b . 

In many experimental situations, the ions move in a plasma with random orientations of d with respect to the 
direction of the polarization of laser field a. The stopping power appropriate to this situation may be obtained by 
carrying out a spherical average over of £(7, a, i?) in Eqs. ( TToT ) and ( ITTb . We find 



5 av (7)«) = s b 



h(7,a) + 3 £2(7, a) 



= S , s ^ av (7, a) 



(21) 



The study of the effect of a radiation field on the SP is easier in the case of low-intensities Wl when a < Ad- 
Then considering in Eqs. (IT&b — (fl~8T> only the quadratic terms with respect to a for the SP S(^, a, d) we obtain 



5(7,0,1?) = Si 



1 "5^( 2C0S 



0+l)Z?( 7 ,O 



(22) 



where 



dx 



(x 2 



1 



Im 



(1 + x 2 w{xh)Y 



(23) 



Taking into account that 7 < 1 and £ S> 1 from Eqs. d22l and (T23l we finally obtain D(7, £) ~ 3/47 2 . It is seen 
that at low-velocities the SP 5(7, a, 1?) decreases with the intensity of radiation field. 

In Figure[T]the quantities 2(7, a, and S av (7, a) are shown vs the intensity parameter a/ Ad of the laser field 
for three values of angles 1? = 0, 1? = 7r/4, ■& = tt/2 and for ujq — 1.2lu p . It is convenient to represent the intensity 
parameter a/Ao in the form o/Ad = 0.18Aq-\/«oW / l/T 1 , where the wavelength (Ao) and the intensity (Wl) of 
the laser field and the density (no) and the temperature (T) of plasma are measured in units /im, 10 15 W/cm 2 , 
10 20 cm -3 and keV, respectively. As an example consider the case when the electron quiver amplitude reaches 
the Debye screening length, a = Ad- For the values of the RF and plasma parameters with Ao = 0.5 /im, no = 
10 18 cm~ 3 , T = 0.1 keV, the above condition is fulfilled at the radiation field intensity Wl = 4.94 x 10 18 W/cm 2 . 

From Figure [T]it is seen that the intense laser field may strongly reduce the SP of the low-velocity ion. And as 
expected the effect of the radiation field is maximal for $ — 0. Note that in this case and at a = Ad the radiation 
field reduces the energy loss Sb approximately by 15 %. For explanation of the obtained result let us consider a 
simple physical model. The stopping power of the ion is defined as S = —(l/v)(dW/dt), where (dW/dt) is the 
averaged (with respect to the period of the radiation field) energy loss rate. We assume that the frequency of the 
radiation field is larger than the effective frequency of the pairwise Coulomb collisions v e s. Also assuming 
that in the low-velocity limit the energy loss of the ion on the collective plasma excitations is negligible and is 
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mainly determined by the Coulomb collisions we obtain (dW/dt) ~ v c ftW . On the other hand u e g ~ V^cff' 
where t> c ff is the averaged relative velocity of the colliding particles. At v < w t h and for vanishing radiation 
field u c ff ~ v t h- However, in the presence of the radiation field the averaged relative velocity of the collisions is 
v c g ~ (v 2 h + v%y/ 2 and increases with the intensity of the laser field. Thus the effective collision frequency v e s 
and hence the stopping power of the ion are reduced with increasing intensity of the radiation field. 

At the end of this section we consider a practical example. Let us consider the stopping of the a-particles in 
the corona of the laser plasma. Although the thermonuclear reactions mainly occur far below the critical surface 
the stopping length of the a-particles is larger than the characteristic length scale of plasma inhomogenei ty and 



some part of the a-particles transfer the energy to the plasma corona before they reach the critical surface (Max 



1982). In the vicinity of the plasma critical density the intensity of the radiation field is very large and the stopping 



capacity of the plasma may be strongly reduced. In this example the typical temperature is T = 10 keV and 
therefore i; Q /i>th = 0.22 (E a = M a V 2 /2 = 3.5 MeV, where E a , M a , v a are the energy, the mass and the 
velocity of the a-particles). For Ao = 0.5 /Ltm, Wl — 2 x 10 17 W/cm 2 , and ujq = uj p ^/2 (the plasma density is 
no = n c /2, where n c is the plasma critical density) we find a ~ Ad- In this parameter regime the radiation field 
reduces the SP of the a-particles by 20 %. 



4 Energy loss of fast ions 



In this section we consider the energy loss of a fast heavy ion moving in a fully degenerate plasma (which means 
that the partially degenerate case could be postponed to a further presentation) in the presence of a radiation 
field. The l ongitudinal dielectric function of t he degenerated electron gas is determined by Lindhard's expression 
( Lindhardi Il954t iLindhard & Wintrier! Il964 ). However, here we consider the simplest model of the dielectric 



function of a jel lium. Previously a plasmon-pole approximation to e(k, to) for an electron gas was used for calcu- 
lation of the SP (Ba sbas & RitchieLfl98^lDeutschLl l995l:lNe"rsisvan & DasLl2000h . In order to get easily obtainable 



analytical results, Basbas & Ritchie ( 1982) employed a simplified form that exhibits collective and single-particle 
effects 

-l ,M 



Im- 



7rajp- — - [S (w 2 - w 2 ) H (k c - k) + 5 (w 2 - uj 2 ) H(k - fc c )] 



e{k,w) 1 "' A/ !24) 

where H(x) is the Heaviside unit-step function, ojk = hk 2 /2m, k c = (2mwj,/fi) 1 ' 2 , and uj p is the plasma 



frequency. The cutoff parameter fc c is determined by equating the arguments of the two delta-functions in Eq 
at k = k c . The first term in Eq. (l24t describes the response due to nondispersive plasmon excitation in the region 
k < k c , while the second term describes free-electron recoil in the range k > k c (single-particl e excitations). Note 
that this approximate dielectric function satis fies at arbitrary k the usual frequency sum rule (IBasbas & Ritchie . 
ll982tlDeutschlll995l:lNersisvan & Dasll2000l) . 

In contrast to the previous section we consider here the fast projectile ion with v > v c (where v c = (jj v /k c 



(fiw p /2m) 1 / 2 ) which justifies the approximation (124-b valid only in this specific case jBasbas & RitchieUl982l) . 

It is constructive to consider first the case of a weak radiation field (k c a < 1) at arbitrary angle $ between v 
and a. In this case the SP is determined by Eq. (Tty. where the field-fr ee SP Sb in the high-velocity limit is given 
by dLindhardl. Il954t ILindhard & Wintheii Il964t beutschl Il986l[l995l) 



Sb 



Z 2 e 2 io 2 



■In 



2mv i 



fkdr, 



Inserting Eq. d24l i into (fl~3T > for the stopping power we obtain 

2Z 2 T, f, , , (k c a) 2 



A 2 



In A 



$ l(A,7) + o $ 2(A,7) sin 2 ?? 



(25) 



(26) 



where Eq = e 2 k 2 = 2fraj p /ao, ao is the Bohrradius, $i = $i c +$i s , $2 = &2c+&2s, A = v/v c , 7 = oj p /ojq < 1 
Also 

3 / 1 \ 3 

7 



1 



4lnA+^ + lVln. 7 



V 



7 



*2c(A, 7 ) = -3 



1 + 7 



$ic(A,7) + 



- - 1 ) In 

7 ) I-7 



2 7 2 A 2 



(27) 
(28) 



$i s (A )7 ) = ^ 



Sf) + - (I3i + Si - ai - rii) 
7 



7 3 V/3 



Vi Si 



atiSi 



1 — A 



(29) 
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^2 S (A,7) = 



1 - 



7 A 2 



47 V ai7?i 



7A 2 ctiS 



4 7 3 A 2 
/8H7i^ + 1 



TJl - S 1 

4 
1 



1 

ai 



- 1 



9 

7A 2 
1 1 
' — + 7" 



-=S I (A 2 + >!) 



(30) 




2 ! = It- w- i - 1 -1 • (3d 

A, 



In Eq. OH n is a positive integer (n = 1, 2, ...)■ The first term in Eq. (l26l i corresponds to the field-free SP (l25T l 
represented in a dimensionless form. The remaining terms proportional to the intensity of the radiation field (a 2 ), 
describe the collective (proportional to $i c; 2c(A, 7)) and single-particle (proportional to $i s; 2s(A, 7)) excitations. 
It should be noted that the stopping power Eq. d26l i is not vanishing only at high- velocities when A ^ 2/ .^7. 

Consider next the angular distribution of the SP at low-intensities of the RF. An analysis of the quantity 
P = (S — Sb)/Sb (the relative deviation of S from Sb) for the proton projectile shows that at moderate velocities 
(A > 2/y/Y) the angular distribution of P has a quadrupole nature. At ^ # ^ ??o(A, 7), where i?o(A, 7) is some 
value of the angle %% the excitation of the waves with the frequencies ojq ± u> p leads to the additional energy loss. 
At $o(A, 7) ^ $ ^ ir/2 the proton energy loss changes sign and the total energy loss decreases. When the proton 
moves at angles 1? = i?o(A, 7) with respect to the polarization vector a the radiation field has no any influence on 
the SP. However, at very large velocities (A 3> 2/^/7) the relative deviation P is negative for arbitrary d and the 
radiation field systematically reduces the energy loss of the proton. 

Let us now investigate the influence of the intense radiation field on the stopping process when v is parallel to 
a. It is expected that the effect of the RF is maximal in this case. From Eqs. ( fTOb and (l24l we obtain 



s = * + £ S £ lE(H J » ( ^' )k ;77Vl 



,n=l 



where A = k c a, P n {x) = (1/ A) (71/7 + x), Q n (x) = (l/A)(n/7 - x), p n = P„(l), q n = Q n (l), and 



(32) 



S = 



A 2 



In A 



dx 



l/A 



Jo W 



is the SP without emission or absorption of the photons. Also we have introduced the notations 

n± = int f kcV^uA = . nt ^ (A ± 1}] ^ 

N = int — — = int 1 — 



(33) 



(34) 



where int(x) is the integer part of x. The quantities a n (A), /3 n (A), <5 ra (A), ?7n(A) in Eq. (1321 are determined by 
Eq. (I3TI 1. We note that in Eq. d32l the terms involving n± and N photons are not vanishing at A ^ I/7 =p 1 and 
A 2/y^7, respectively. Similarly the SP (l33l is not vanishing at A ^ 1. 

The first term in Eq. (l33~t describes the collective excitations while the second term corresponds to the single- 
particle excitations. From Eq. (l33l it is seen that So oscillates with the intensity of the laser field. However, the 
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ka ka 

c c 



Figure 2: (Left panel) the ratio R(a) — S(o)/Sb as a function of dimensionless quantity k c a at v = 8.6tv, ujq = 1.2ui p 
(solid line), ujq = l.&u p (dashed line), uio = 2ui p (dotted line), uio = 3uu p (dash-dotted line). Thin solid line corresponds to 
Ro(a) = So(cl)/Sb (see Eq. | |33I >). (Right panel) same as in left panel but at ujq — 1.2ui p , v — 3v c (solid line), v — 7v c 
(dashed line), v — llv c (dotted line), v — 17v c (dash-dotted line). 



radiation field suppresses the excitation of the collective and the single-particle modes and the SP So is less than 
the field-free SP Sb- As follows from Eq. (l33l at high-intensities of the RF the SP So is close to zero when 
A/\ ~ fi m (or alternatively at j(ve/v) ~ /i m ) with m = 1, 2, . . ., where fj, m are the zeros of the Bessel function 
Jo(Hm) — (/ii = 2.4, H2 = 5.52, /i3 = 8.63. . . ). Then the energy loss of the ion is mainly determined by the 
other terms in Eq. (l32l i and is stipulated by excitation of plasma waves with frequencies tiluq ± oj p . The first and the 
last pairs of terms in Eq. (|32| | describe the excitation of the collective and single-particle modes, respectively, with 
emission or absorption several photons. The number of photons (n±, N) involved in the process of the inelastic 
interaction are determined by the energy-momentum conservations (see the arguments of the delta-functions in 
the dielectric function (l24l ). 

The results of the numerical evaluation of the SP (Eqs. (l32l > and d33"T >) are shown in Figure [2] where the ratio 
R(a) — S(o,)/Sb is plotted as a function of the laser field intensity (k c a = 5.38W^ 2 uiQ 2 r s 3 ^ 4 , where r s is the 
Wigner-Seitz density parameter and Wl and ujq are measured in units 10 15 W/cm 2 and 10 16 sec -1 , respectively). 
For instance, for Al target with 7- s = 2.07, haj p — 15.5 eV, and v c = 1.2 x 10 s cm/sec. From Figure |2]it is seen 
that the SP exceeds the field-free SP and may change sign due to plasma irradiation by i ntense (k c a ^> 1) laser 
field. Similar properties of the SP has been obtained previously for a classical plasma dNersisvan & Akopvan 
1999). However, due to the higher density of the degenerate electrons (in metals typically uq ~ 10 23 cm 3 ) 



the acceleration rate of the projectile particle is larger than similar rate in the case of a classical plasma. The 
acceleration effect occurs at ve/v ~ Mm/7 (with m = 1,2,...) when the SP So nearly vanishes. It should 
be noted that in the laser irradiated plasma a parametric instability is expected dSihnl 1 1973b with an increment 
increasing with the intensity of the radiation field. This restricts the possible acceleration time with stronger 
condition than in the case of a classical plasma. Finally, let us note that the effect of the enhancement of the SP 
of an ion moving in a laser irradiated plasma is intensified at smaller frequency (Fig. [2] left panel) of the radiation 
field (ojo — 0J p but uiq > u P ) or at larger incident kinetic energy of the projectile ion (Fig. [2 right panel) when the 
numbers n± and N of the photons involved in the inelastic interaction process are strongly increased (Eq. (l3~4l). 



5 Summary 

In this paper, within RPA we have investigated the energy loss of a heavy point-like ion moving in a laser irradiated 
plasma. In the course of this study, we derived a general expression for the SP which has been also simplified in 
the limit of weak RF. As in the field-free case, the SP in a laser irradiated plasma is completely determined by the 
dielectric function of the plasma. We have considered two somewhat distinct cases of the slow- and high-velocity 
ion moving in a classical and fully degenerate electron plasma, respectively. At low-velocities the RF leads to 
the strong decrease of the energy loss. Physically, this is due to the strong reduction of the effective frequency 
of the pairwise Coulomb collisions between projectile ion and the plasma electrons. At high velocities the RF 
may strongly increase the SP. This effect is more pronounced when the laser frequency approaches the plasma 



frequency in agreement with PIC simulations dHu et al. . 2011 ). Moreover, at high-velocities and in the presence 



of the intense RF an ion projectile energy gain is expected when the quiver velocity of the plasma electrons exceeds 
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the ion velocity. The analysis presented above can in principle be extended to the case of a partially degenerate 
plasma as well as to the case of light ion projectiles and also electrons and positrons when the effect of the intense 
RF on the ion cannot be neglected anymore. We intend to address these issues in our forthcoming investigations. 
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